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Abstract 

A reusable launch vehicle control problem during ascent is addressed via multiple-time 
scaled continuous sliding mode control. The proposed sliding mode controller utilizes a two-loop 
structure and provides robust, de-coupled tracking of both orientation angle command profiles 
and angular rate command profiles in the presence of bounded external disturbances and plant 
uncertainties. Sliding mode control causes the angular rate and orientation angle tracking error 
dynamics to be constrained to linear, de-coupled, homogeneous, and vector valued differential 
equations with desired eigenvalues placement. Overall stability of a two-loop control system is 
addressed. An optimal control allocation algorithm is designed that allocates torque commands 
into end-effector deflection commands, which are executed by the actuators. The dual-time scale 
sliding mode controller was designed for the X-33 technology demonstration sub-orbital launch 
vehicle in the launch mode. Simulation results show that the designed controller provides robust, 
accurate, de-coupled tracking of the orientation angle command profiles in presence of external 
disturbances and vehicle inertia uncertainties. This is a significant advancement in performance 
over that achieved with linear, gain scheduled control systems currently being used for launch 


vehicles. 



Introduction 


Flight control of both current and future reusable launch vehicles (RLV) in ascent mode 
involves attitude maneuvering through a wide range of flight conditions, wind disturbances, and 
plant uncertainties. A flight control algorithm that is robust to changing flight conditions, 
disturbances, and vehicle uncertainties would be an improvement over current RVL flight 
control technology. Sliding Mode Controller (SMC) is an attractive robust control algorithm for 
RLV ascent designs because of its inherent insensitivity and robustness to plant uncertainties and 
external disturbances 1 ' 3 . 

A SMC design consists of two major steps: (1) a sliding surface is designed such that the 
system motion on this surface exhibits desired behavior in the presence of plant uncertainties and 
disturbances; and (2) a control function is designed that causes the system state to reach the 
sliding surface in finite time and guarantees system motion in this surface thereafter. The 
system’s motion on the sliding surface is called sliding mode. Strict enforcement of the sliding 
mode typically leads to discontinuous control functions and possible control chattering effects 4 . 
Control chattering is typically an unwanted effect that is easily eliminated by continuous 
approximations of discontinuous control functions 4 ' 3 , or by continuous SMC designs 6 . 

An example of the application of SMC is shown in the spacecraft attitude control work 
performed by Dwyer et al 7 . This work utilizes the natural time-scale separation, which exists in 

g 

the system dynamics of many aerospace control problems . 

The RLV SMC architecture, which is developed in this paper, is a two loop structure 0 
that incorporates backstepping techniques 10 . In the outer loop, the kinematics equation of angular 
motion is used with the outer loop SMC to generate the angular rate profiles as virtual control 
inputs to the inner loop. In the inner loop, a suitable inner loop SMC is designed so that the 


2 



commanded angular rate profiles are tracked. The inner loop SMC produces roll, pitch and yaw 
torque commands, which are optimally allocated into end-effector deflection commands. 
Multiple time scaling (multiple-scale) is defined as the time -constant separation between the two 
loops. That is, the inner loop compensated dynamics is designed to be faster than the outer loop 
dynamics. The resulting multiple-scale two-loop SMC, with optimal torque allocation, causes the 
angular rate and the Euler angle tracking errors to be constrained to linear de-coupled 
homogeneous vector valued differential equations with desired eigenvalues placement. 

Simulation results are presented that demonstrate the SMC’s effectiveness in causing the 
X-33 technology demonstration Launch vehicle, operating in ascent mode, to robustly follow the 
desired profiles. The resulting design controls large attitude maneuvers through a wide range of 
flight conditions, provides highly accurate tracking of guidance trajectories, and exhibits 
robustness to external disturbances and parametric uncertainties. 

Equations of RLV dynamics 

The dynamic equations of rotational motion of the rigid-body RLV are given by the Euler 
equation in the body frame 

(j 0 + Aj)cu = — Q(j 0 + Aj)w + T + d , (1) 

where J 0 6 R 3x3 is the nominal inertia matrix, AJ e R 3 * 3 is an uncertain part of the inertia 
matrix, caused by fuel consumption and variations of particular payloads from a nominal one, 
o> = [ p q r] T is the angular rate vector, T = {L,M,iV} r is the control torque vector, 
d = {z.„Af„./U r is the external disturbance torque vector. The matrix Q, is given by: 
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a), 0 -ft), . (~ ) 
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The RLV orientation dynamics are described by the kinematics equation on the body 


y = R(y)£0, 
where 



1 

tan 6 sin cp 

tan 9 cos cp 

R 00 = 

0 

cos cp 

- sin cp 


0 

sin cp 
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cos 9 

cos 9 

y = [cp 

Q 

vf 



The RLV equations of the translational motion are given in the body frame as 

v = Hco - g<!> + — F (5) 

M 

where 

0 —v. v v 

H = v z 0 v t , ( 6 ) 

_-v v v ( 0 _ 

v = [v t , v v , v. is the linear velocity vector, M is the mass of RLV, F is the vector of external 

forces acting on RLV, <t> = [sin cos0sin<p cos#cos<p] r , and g= 9. Sim /s' is a 

gravitational constant. The total mass expenditure onboard the RLV in the launch mode is 
described as follows: 


M-m dry .+m f , m f =m JO -m f t, m { - 


m 0 , if m f >0 
0, ifm f =0 
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where m f and m f0 are the current and initial fuel masses, respectively, m dry is a mass of the 
RLV without fuel (dry), and rhj is a rate of fuel consumption. 

The control torque T is generated by the aerodynamic surfaces and rocket engines. This 
is described by the equation 
T = D(.)<5 , 

where D(.)e R 3 '" is a sensitivity matrix calculated on the basis of table lookup data, <5e R'' is 
the vector of aerodynamic surface deflections and differential throttles of the rocket engines. 

Electromechanical actuators deflect the aerodynamic surfaces (aerosurfaces), and thrust- 
vector-control valves that throttle the rocket engines. The actuator dynamics are assumed to be 
much faster than the dynamics of the compensated RLV flight control system. 

Problem Formulation 

The general problem formulation is to determine the actuator deflection commands <5 ( 
such that the commanded Euler angle profiles y c = [<p c 0 C i/L ] ? are r °bustly followed in the 
presence of bounded disturbance torques d and plant inertia variations AJ . This problem is 
addressed as follows: 

• A control law is designed in terms of the control torque command vector T f , and 

• An optimal control allocation matrix B, is designed to map the control torque command 
vector T to actuator deflection commands. That is, 

S c = B/T c . ( 9 ) 

• The “fast” actuator dynamics are neglected at the stage of the controller design. They must be 
used during simulations to validate the designed controller. 
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Specifically, the control problem for the RLV in launch mode is to determine the control 
torque input command vector T in the given state variable equations 

'(J 0 + Aj)w = -0.(3 0 + Aj)w + T, + d, 

<y = R(y)G>, (10) 

y = Y 

such that the output vector y asymptotically tracks a command orientation angle profile 
y r = Y c ■ That is, 

lim|y lf (/)-y i (/)|| = 0 Vi = l,3. (11) 

t — >°° 

The desired RLV performance criteria is to robustly track both the commanded Euler 
angle guidance profiles y c and the real-time generated angular rate profiles (i ) r , such that the 
motion for each quantity is described by a linear, de -coupled, and homogeneous vector valued 
differential equations with given eigenvalues placement. 

Sliding Mode Controller Design 

Tracking of the commanded Euler angles guidance profile y, (^ e primary 

objective) is achieved through a two loop SMC structure 9 . The cascade structure of the system in 
eq. (10) and the inherent two time scale nature of the RLV flight control problem are exploited 
for design of a two loop flight control system using continuous SMCs in the inner and outer 
loops. The outer loop SMC provides angular rate commands 0) c to the inner loop. The inner loop 

SMC provides robust de-coupled tracking of angular rates (t ) . Together, the inner and outer loop 
SMCs form a two loop flight control system that achieves de-coupled asymptotic tracking of the 
command angle guidance profile y c . Similar structure was developed" for an aircraft control 

system using a dynamic inversion algorithm. 
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Outer Loop Sliding Mode Controller Design 

The outer loop SMC takes the vehicle's angular rate vector (0 as a virtual control £0 and 

uses the kinematics eq. (3) to compensate the Euler angle tracking dynamics. The motion of the 
outer loop compensated tracking error dynamics, with desired bandwidth, is constrained (by 
proper control action) to a sliding surface of the form 

a = y t +K 1 jY t .£/r , creR\ (12) 

0 

where y e =Y C ~Y and K, = diag{k u }, K, e R 3x3 . The gain matrix K t is chosen so that the 
output tracking error y e exhibits a desired linear asymptotic behavior on the sliding surface 
(<T=0). 

The objective of the outer loop SMC is to generate the commanded angular rate vector 
£U ( . (which is passed to the inner loop SMC) necessary to cause the vehicle’s trajectory to track 

the commanded Euler angles guidance profile y r . In other words, the virtual control law 0) c is 
designed to provide asymptotic convergence (with finite reaching time) of the system’s eq. (3) 
trajectory to the sliding surface a = 0 . Dynamics of the sliding surface in eq. (12) are described 
as 

d = y t . -R(y)(o r +K,y r . ( 13 ) 

The outer loop SMC design is initiated by choosing a candidate Lyapunov function of the 

form 

V =—a T a >0, (14) 

2 

whose derivative is shown as 


V = a T a =<J r [y c -R(y)w c +K,yJ. 


(15) 



To ensure asymptotic stability of the origin of the system in eq. (13), the following 
derivative inequality of the candidate Lyapunov function is enforced 1 4 

3 

V =-po T SlGN{o) = -pY J \o,[ P> 0. ( 16 ) 

1=1 

Considering equality (16), the required angular rate command (O c to ensure asymptotic 
stability is defined as 

(d c =R-'(y)[y c +Kj,] + R-'(y)PSIGN(<j), (18) 


where SIGN (a) = [s/gn(<7, ), sign(o 2 ), sign(o 3 )f , and p > 0 . 

The sliding surface shown in eq. (12) will reach zero in a finite time 1 4 defined as 


K 


= max 

*6(1.31 


K <o)| 

J 

P 


(19) 


where p > 0 , and t r is the design parameter describing the sliding surface reaching time. The 


value of p is calculated using inequality (19). 

The angular rate command profile in eq. (18) is discontinuous; as a result, the control 
actuation will chatter during system operation on the sliding surface eq. (12). Mechanically and 
electrically, this chattering is an unwanted effect. Moreover, a discontinuous profile cannot be 
accurately tracked in the inner loop. To solve this problem, the discontinuous term ( SIGN(cr)) 
in eq. (18) is replaced by the high-gain linear term 

f -\T r s 


K n o = 


<7, cr ; <7 


*K 0 = diag\ — \ Ve, >0 VI, 3, 


( 20 ) 




where e, > 0 VI, 3 define the slopes of the linearized function 4,5 . 

After substitution of eqs. (20) and (18) into eq. (13), the sliding surface dynamics are of 
the form 
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6 - -pK 0 cr . 


( 21 ) 


Eq. (21) is globally asymptotically stable, since p >0 and K 0 is positive definite, and 
the equilibrium < 7=0 will be reached asymptotically. Moreover, in a close vicinity of < 7 , = 0 , 
Vi = 1,3, the tracking error y t will exhibit de-coupled motion in accordance with eq. (14). 

Inner Loop Sliding Mode Controller Design 

The purpose of the inner loop SMC is to generate the vehicle torque command vector T 

necessary to track the given commanded angular rate profile (O c . In addition to solving the inner 
loop tracking problem defined as 

lim||ft>, r (/) - d) i (r)|| =0, Vi = 1,3 , (22) 

SMC causes the system to exhibit linear decoupled motion in sliding mode. 

The motion of the inner loop compensated dynamics, with desired bandwidth, is 
constrained (by proper control action) to a sliding surface of the form 

t 

s = w e +K 2 j(o e dr, se R 3 , ( 23 ) 

0 

where co e =(O r -w and K 2 = diag{k 2 }, K,eR 3x3 . The angular rate tracking error (o e exhibits 

a desired linear asymptotic behavior while operating on the sliding surface (i = 0) and with 
proper choice of K, . The inner loop sliding mode dynamics in eq. (23) are designed faster than 
the outer-loop sliding mode dynamics in eq. (14) in order to preserve sufficient time scale 
separation between the loops. 
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The command torque T r is designed to provide asymptotic convergence (with finite 

reaching time) of the system's eq. (1) trajectory to the sliding surface s- = 0 . Dynamics of the 
sliding surface eq. (23) are described as follows: 

j = d> 4 .+(j 0 +Aj)“ , Q(j 0 +Aj)ft)-(Jo+Aj)‘ , T f -(j 0 +Aj)" l d + K 2 fl),. (24) 

d\\A3\\ 

Considering that the inertia matrix J 0 + AJ is positive definite and assuming " negligibly 
" at 

small, the inner loop SMC design is initiated by choosing the candidate Lyapunov function of the 

form 


V = — s 7 (j 0 + Aj)s > 0. 


(25) 


whose derivative is shown as 


V = s T 


-Ajs + (j 0 + Aj)w r +(J 0 + Aj)K 2 cu f +Q(j 0 + AJ)m-T. -d 


(26) 


To ensure asymptotic stability of the origin of the system eq. (24), the following 


derivative of the Lyapunov function candidate is enforced 

3 

V <-r]s T SIGN (s) = t \, V > 0. 


1-4 


(27) 


Further, the sliding surface eq. (23) will reach zero in a finite time 1-4 defined by 

s,(0)| 


r < max 


(28) 


ie[U] fj 


where t) > 0 , and t r is the design parameter describing the sliding surface reaching time. The 
value of rj is calculated using inequality (28). 

Considering inequality (27), the required torque command T c to ensure asymptotic 
stability is defined as 
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T =J 0 o) ( . +J 0 K 2 o), +nj 0 (0 + ps + pSlGN(s). (29) 

One should note that the SMC eq. (29) is independent of uncertainties AJ, AJ and disturbances 
d . Substituting eq. (29) into eq. (26) yields 

U = +-^ r Aj^+5 7 '(AJd) ( . +AJK : £O t , + £2AJft>-d) . (30) 

i=i 2 

Using the inequality 4 s r AJs < A||s|| : , where A is greater or equal to the maximum eigenvalue of 
the matrix AJ , and assuming 

|(Ajft) ( .) i | < a n |(AJK 2 G>,) ( .| < b { , |(aAj£u),|<c,, |(d),|<L, Vi = U, (3D 

then inequality (30) is rewritten as follows: 

V<-£(p-a, -b, -c, -L ( )|5,.|-(p-A)|5||'. (32) 

i=\ 

Inequality (32) is rewritten to enforce inequality (27) after selecting fx > A . This is, 

V < -£ (p - a, - 6, - c, - L, )|s, | < \s, | . (33) 

i=I i=l 

The value of parameter p is identified in accordance with inequality (33) as follows: 
p>a t +b t + c, + L, + T) (34) 

Remark. The constants a l ,b l ,c i ,L i Vi = 1,3, are derived by estimating the upper bounds of 

inequalities (31) within a reasonable flight domain. 

To avoid chattering, the SMC eq. (29) is implemented in a continuous form 
T c =J 0 (b c +J 0 K 2 w e +QJ 0 a) + ixs + pSAT(rs), (35) 

where 
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r = ding < 


C 


V/ = 1,3, SAT(Fs) = \ sat sat— -,sat - 3 


_fl_ 

c, 


S 2 

cT 


^3 

C 3 


•v, 

A-flf-T-' 
V / 


i- a > c, 

s, 


if s t < C, (36) 


c 

- 1, 'A < -C< 


Two-loop stability analysis 

The outer loop continuous SMC in eqs. (18) and (20) provides global asymptotic stability 
of the equilibrium point a = 0 in accordance with eq. (21). However, the virtual control signal 
(o c defined by eqs. (18) and (20) is tracked in the inner loop with a tracking error (i) e . Using 
backstepping techniques 10 , the compensated dynamics of the outer loop sliding surface in eq. 
(21) must be rewritten as follows 

6 = -pK 0 a + R(y)co e . ( 37 ) 

Initially, suppose the virtual control signal (O c is tracked in the inner loop by the 
discontinuous SMC in eqs. (29) and (34). This SMC provides global asymptotic stability of the 
equilibrium point 5 = 0, since the negative definiteness of the Lyapunov function derivative in 
eq. (33) is maintained (given sufficient control power) in the presence of inertia matrix 
uncertainties and external disturbances. Furthermore, the equilibrium point 5=0 is reached in 
the finite time t - T r in eq. (28). So, Vr > z r the inner-outer loop compensated dynamics are 
described as follows: 


1 

q. 

1 

— 

X o 

1 

1 

l ' 

?! * 

o 

-K(.)K, 

-/?(.) k 2 

^ r 5 ' 

co e =-K 2 z. 

(38) 

J - 


0 

0 

K 2 

Z 




In accordance with eq. (38) and that matrix K, is chosen positive definite, the 
equilibrium z - 0 (and = 0 ) is globally asymptotically stable. The equilibrium a = 0 is also 
globally asymptotically stable, since the matrix K 0 is chosen positive definite and p > 0 . 
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Furthermore, the equilibrium y e = 0 is globally asymptotically stable, since the matrix K, is 
chosen positive definite. Consequently, the system in eq. (38) has a globally asymptotically 
stable equilibrium {y e ,o,zY = 0 . 

The y t ar *d error dynamics are de-coupled using the multiple time-scale concept. 
This is accomplished by selecting positive definite matrices K, and K-, such that 
min{Re[e/ < g(K,)]}» max{Re[<?/g(K, )]}, and defining the elements of the diagonal positive 

definite matrix K 0 = dj'flgj— | V/ = 1,3 and the scalar p> 0 such that 

min{Re[<?/g(pK 0 )I}» max{Re[etg(K, )]}. This approach causes the z- and a - transients to 
die out first followed by the lim y e - 0 according to y e - -K . 

Now, suppose the virtual control a> c is tracked in the inner loop by the continuous SMC 
in eqs. (35) and (36). The derivative of the Lyapunov function in eq. (26) is calculated as 
follows: 


V< 


‘f N>£ 

t =i 

3 f i 


-» 


pY 

C 


(39) 


~ a > ~ b i ~ c i 


if $ 


The boundary layer ^1 = ^, V/ = l,3 will be reached in finite time t = T r , since V 7 <— 

1=1 


V s t >^ 4 ' 5 . The resulting inner-outer loop compensated dynamics, Vr > x r , are described as 


follows: 



■yl f-K,-pK 0 -R(.)K 2 K(.)Ty e l ro 
a _ 0 - pK 0 -/?(.) K : /?(.) cr 0 

Z “ 0 0 -K 2 1 3 x 3 2 + 0 

S — A 41 ~ A 4 -> — A 43 ~ ^44 _ _ _^-»5 

where 

A 41 = [(l-(I + E)-')?' , (r)(K l + /?(.)/?- , (/))-(J 0 +AJ)- l QAJ«- 1 (y>>C 1 , 

A 42 = [p(l - (I + E)- 1 )r~' (y)(K, + pK 0 + R(y)R~ l (/))- (J 0 + AJE'^AJ/T 1 (y)jK 0 , 

A 43 = [(i - (I + E)- 1 >?-' (y )(K, + pK 0 )R(y ) + (i - (I + E)~%_ - (J 0 + AJ)' 1 ^AJ Jk 2 , 

A„ = (J 0 + AJ)- 1 (pi + pH + (Jo + AJ)' 1 £>AJ-(l-(I + E)- 1 >?-' (y)(K, + pK 0 )/?(y) -(i-(I + E)~')k 2 , 
A 45 = -(J 0 +AJ)"‘, A 46 = (l-(I + E)-'>?- 1 (y), 

A 47 = -((i - (I + E)' 1 ]/T‘(y)tf(y)- (J 0 + AJ)' 1 ^Aj)/?' 1 (y) , E = J 0 -AJ . 

The first three equations in eqs. (40) are the same as eqs. (38), with the exception of the 
coupling of the bounded signal 5 . The system in eq. (38) has a globally asymptotically stable 
equilibrium {y t , a, zY = 0 . Hence, the forced responses described by the first three equations in 

eqs. (40) are bounded with boundary widths proportional to £, Vi = 1,3 . 

In particular, if AJ=0, then A 41 = A 42 = A 43 = A^ = A 47 = 0 , A 43 =J 0 and 

A w =J 0 '(pI + pH • Matrix A 44 is obviously positive definite and the equilibrium 5 = 0 is 
globally asymptotically stable. Further, if limd(r) = C = const , then 

lims = (pi + pT)’ 1 C = const and limz = K 2 ‘(pl + pr)~‘C = const . It is easy to show that in 

/— » oo l — ><» 

this case limcr = limy^ = 0 . 
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The following are conclusions derived from the two-loop stability analysis of the 
continuous inner and outer loop SMCs operating in the presence of bounded inertia AJ and 
external torque d disturbances: 

• The asymptotic transient behavior of the inner and outer loop tracking errors y e an< ^ > 

which are controlled by the continuous outer loop SMC and the discontinuous inner loop 
SMC, are de-coupled in presence of bounded AJ and external disturbance torque d . 

.• The tracking errors y t and 0 ) e , which are controlled by the continuous outer loop SMC and 

the continuous inner loop SMC, are bounded in presence of bounded AJ and external 
disturbance torque d . 

• The asymptotic transient behavior of the inner and outer loop tracking errors y t and co e , 
which are controlled by the continuous outer loop SMC and the continuous inner loop SMC, 
are de-coupled in presence of external disturbance torque d: limd(r) = C = const and 

completely known matrix of inertia J . 


Control Allocation 

The command torque T is related to the deflection command vector 8 C e R" as follows: 

T = D(.)<5 ( . ( 41 ) 

Assuming rankD(.) = rank\Q(.), T r ]= 3, it is obvious that given a T. profile, there exist an 
infinite number of corresponding S ( . profiles that satisfy eq. (41). This lends itself to the 


following optimization problem 

min f (8 C ) = 8*Q8 C , fl: {(5 C |D(.)5 C -T c =o}, 

which is minimized by the solution 


( 42 ) 
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(43) 


8, =Q-‘D(.) t [d(.)Q' 1 D(.) t }' , T c . 

Eq. (43) is re-written as 

S. = B.T , 

where 

B v = Q~'D(.) t [d(.)Q~ 1 D(.) t ]” ( 45) 

is the control allocation (gain-mixing) matrix, Q is the weighting matrix, and D(.) is the torque- 
to-deflection partial derivative matrix or sensitivity matrix. 

Design and Simulation of the X-33 RLV Control System in the Launch Mode 

The X-33 technology demonstration RLV controller design is considered in the launch 
mode. The 6DOF high fidelity mathematical model includes the control deflection vector 
5 6 R 11 that consists of the following components: 8 i and S 2 are deflections of the right and left 
flaps (in radians); S 2 and <5 4 are deflections of the right and left inward elevons (in radians), 
and S 6 are deflections of the right and left outward elevons (in radians); 5 7 and 8 S are 
deflections of the right and left rudders (in radians); <5 9 , 5 10 , and <5 U are pitch, roll and yaw 
differential throttles (in % of power level) of the aerospike rocket engine. Flight guidance tables 
for the Euler angle command profiles y c are used. High fidelity sensor mathematical models 
with 0.02 s time delay as well as realistic table look up wind gust profiles that create a 
disturbance torque are also used. 

The command profiles co c and T c are generated by the outer and inner loop SMCs in 
eqs. (18), (20) and (35). The deflection command vector 8 C is calculated in eq. (9) and executed 
by the actuators, whose dynamics are contained in high fidelity models used in the simulations. 
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In the calculation of the control allocation matrix B„ in eq. (45), the Q 1 matrix is chosen to be 
diagonal. The diagonal terms are equal to square maximum deflections 8'^ , which provide 
both a dimensionless performance index f (5 r ) = SjQ8 c and a uniform penalty function for each 
command 5. The control torque T is calculated by eq. (8) and used in eq. (1) during the 
simulation. The time-varying sensitivity matrix D(.) is obtained by linearizing through the 
aerodynamic and propulsion system data bases at various flight conditions. 

The outer loop continuous SMC is designed in accordance with eqs. (12), (18) and (20). 
The following outer loop SMC parameters are chosen: 

'1.0 0 Ol [0.4 0 0" 

p = l, K 0 = 0 0.8 0 , K, = 0 0.4 0 , (46) 

0 0 1.5 [° 0 °- 4 

where the matrix K, provides a 0.4 rad Is given bandwidth for the compensated outer loop. The 

inner loop continuous SMC is designed in accordance with eqs. (23), (34) and (35). The 

following inner loop SMC parameters are chosen: 

"3.8 -10 6 0 0 1 [1.88 0 O' 

p = 0 3.8 TO 6 0 N-m, K 2 = 0 1.88 0 , (47) 

0 0 5.0 10 6 .o 0 1.88 

where the matrix K, provides a 1.88 rad /sec given bandwidth for the compensated inner loop, 
and p = 0 is chosen since the matrix AJ 0 is a negative definite. 

The 6DOF time simulation results are presented in figures 1-18. Figures 1-3 and figure 
7 show overall tracking performance of the Euler angle guidance profiles in the outer loop, 
which is obviously decoupled and very accurate. Figures 4-6 show reasonably accurate tracking 
the angular rates command profiles in the inner loop. Figures 14 and 15 demonstrate torque 
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command profiles T,.. These commands are allocated into actuator deflection commands 5 C and 

executed by the actuators. Figures 8-13 demonstrate corresponding deflections ot the 
aerodynamic surfaces and engine differential throttles that are far from saturation. Figures 16-18 
show that tracking motion happens in the close vicinities of the sliding surfaces. This means that 
the robust accurate real sliding mode exists in the X-33 control system by means of the two loop 
multiple scale SMC in the presence of wind disturbances and inertia uncertainties. 


Conclusions 

Employing a time scaling concept in the inner (angular rates) and outer (orientation 
angles) loops, a new two-loop continuous SMC was designed for the RLV . This SMC provides 
de -coupled performance in the inner and outer loops in presence of bounded external 
disturbances and plant uncertainties. Stability of the two-loop control system is analyzed. A 
control allocation matrix is designed that optimally distributes the roll, pitch, and yaw torque 
commands into end-effector deflection commands. 6 DOF high fidelity simulation show that the 
multiple scale SMC provides robust, accurate, de -coupled tracking of the commanded Euler 
angle guidance profiles for the X-33 RLV in the ascent mode in the presence ot external 
disturbances (wind gusts) and plant uncertainties (changing matrix of inertia). 
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